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ABSTRACT 

We present high spatial resolution spectroscopic measurements of dynamic fibrils (DFs) in the 
Ca II 8662 A line. These data show clear Doppler shifts in the identified DFs, which demonstrates 
that at least a subset of DFs are actual mass motions in the chromosphere. A statistical analysis 
of 26 DFs reveals a strong and statistically significant correlation between the maximal velocity and 
the deceleration. The range of the velocities and the decelerations are substantially lower, about a 
factor two, in our spectroscopic observations compared to the earlier results based on proper motion 
in narrow band images. There are fundamental differences in the different observational methods; 
when DFs are observed spectroscopically the measured Doppler shifts are a result of the atmospheric 
velocity, weighted with the response function to velocity over an extended height. When the proper 
motion of DFs is observed in narrow band images, the movement of the top of the DF is observed. 
This point is sharply defined because of the high contrast between the DF and the surroundings. The 
observational differences between the two methods are examined by several numerical experiments 
using both numerical simulations and a time series of narrow band Ha images. With basis in the 
simulations we conclude that the lower maximal velocity is explained by the low formation height of 
the Ca IR line. We conclude that the present observations support the earlier result that DFs are 
driven by magneto-acoustic shocks exited by convective flows and p-modes. 
Subject headings: Sun: chromosphere — techniques: spectroscopic — Sun: atmospheric motions 



1. INTRODUCTION 

The dynamical nature of the chromosphere is obvious 
when the Sun is imaged in the line center of strong chro- 
mosp heric spectral lines, most commonly in the Ha line 
(e.g. . Ivan Noort fc Rouppe van der Voortll2006f) . One of 
the dominating features is a vast numb er of thin (0.2 - 
1 Mm) omnipresent jet like structures (jBeckersI 119681 ). 
On the quiet solar limb they are commonly known as 
spicules, while on the quiet disk they are often called 
mottles, and finally in active regions they are known 
as active region fibrils or dynamic fibrils (DFs). The 
nomenclature can be confusing, but there are strong 
indicatio ns that these structure s are physically c losely 
related (iTsiropoula et al.l 119941: iSuematsu et al] 119951: 
Christopoulou et al.l 12001 : iRouppe van der Voort et al.l 
2007^^ 

Recently, a combination of high-resolution observa- 
tions and advanced numerical modeling have shown that 
DFs are most likely driven by shocks that form when pho- 
tospheric oscillations leak into the chromosp here along 
inclined flux tubes (Suematsu 1990; De Pontie u et al.l 



[20C)4 HHansteen et al. 2006; Dc Pontieu ct al. 20M)- The 
inclination of the magnetic field lowers the acoustic cut- 
off frequency sufficiently to allow p-modes with the 
dominant low frequencies to propagate along flux tubes 
(lMichalitsanoslfl973HBei fc Lerovl[l97l . 

These insights into the formation of DFs have be- 
come possible because of recent developments in obser- 
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vational techniques, such as bigger telescopes combined 
with real time wavefro nt corrections by adaptive op- 
tics (A O) systems (e.g. iRimmel j l2000t (Scharmer et al. 
2003b), and post-processing methods (e.g. lvon der Lueh" 



1995 



van Noort et al.ll2005( l which have made observa- 



tions of these jet structures much more reliable. 

These developments have spurred several authors 
to focus on the detailed understandin g of DFs (e.g., 
De Pontieu et al.ll2003L 12004,2005. 20071: iTziotziou et all 



2004 iHansteen et al.ll2006ride Wiin fc de Pontieul 120061: 



Koza, et al.l 120071 : iHeggland et al.ll2007t l. One of the im- 
portant results of this work is that the DFs are driven by 
and can channel photospheric oscillations into the chro- 
mosphere and the corona^ 

In two pap e rs Hansteen et al.l (|2006f ) and 
iDe Pontieu et al.l (|2007f l used high spatial and high 
cadence observations of the Ha line center together 
with realistic simulations to investigate the nature of 
DFs. One of their conclusions was that the DFs follow 
parabolic paths along their axis with decelerations lower 
than the solar gravitational deceleration. Previous 
observations did not allow an accurate determination 
of the nature of the trajectory because of lower quality 
data and line-of-sight (LOS) effects that were dif ficult 
t o estimate (|Nis hikawa 1988; Suema tsu et al.iri995l). 

IHansteen eTal.. (2006.) and De Pontieu et al.l (|2007l ) 
further report regional differences between DFs observed 
in two different plage areas. The regional differences are 
explained by different inclination angles of the magnetic 
fields in the two regions. In this way the magnetic topol- 
ogy of the solar atmosphere works as a filter, where only 
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Fig. 1. — Spectrograms of the Ha line (panel a) and the Ca II 8662 A line (panel b). Notice the highly dynamical line center of the 
Ca line. The corresponding MFBD processed slitjaw image (panel c) and the narrow band Ha DOT image (panel d). The position of the 
spectrograph's slit is marked with a line in the DOT image. 



waves with certain periods can leak through. The simu- 
lations, spanning from the upper convection zone to the 
corona, reproduce the observed correlations between the 
maximum velocities and decelerations in DFs leading to 
the conclusion that DFs are formed by chromospheric 
shocks driven by global p-modes and convective flows. 

In this paper we add to the understanding of the DFs, 
by analyzing high spatial and high cadence spectrograms 
of the Ca II 8662 A line, put into context by simultane- 
ous Ha spectrograms and narrow-band images. Further- 
more, the observational results are compared with the 
numerical simulations of iHansteen et all ()2006f ). 

In § [2] we describe the observing program and instru- 
mentation. The data reduction method is described in 
§[3 In §|4]we show the results of the observations. The 
main observational errors are discussed in § [5l To get a 
better understanding of the observations we present sev- 
eral numerical experiments in §[6l Finally, we summarize 
the results in § [71 

2. OBSERVING PROGRAM AND INSTRUMENTATION 

The data presented in this paper were obtained in a 
co-observation ca mpaign between the Sw edish 1-m Solar 
Telescope (SST) (IScharmer et al.ll2003a^ an d the Dutch 
Open Telescope fDOT) (jRutten et al.ll200l . The spec- 
trograms were obtained with the SST, using the TRIP- 
PEL spectrograph^ in combin ation with the ada ptive 
optics system at the SST (.Scharmer et al.l [2003b[ ) . In 
short the TRIPPEL spectrograph is a multi port spec- 
trograph, using an echelle grating with a blaze angle of 
63.43°. The theoretical spectral resolution is 240000. In 
our case we observe in two ports to obtain simultaneous 
Ha and Ca II 8662 A spectra. We observe at —1.06° 
away from the blaze angle in order 34 and 26 respec- 
tively. To remove light from other orders we use two 
standard filters centered at 6562 and 8680 with FWHM 
of 43 and 101 A respectively. These filters are placed 
in front of two Megaplus 1.6 cameras (KAF-0401 Im- 
age sensor). With this setup the spatial pixel size is 

|http://dubshen.astro.su.se/wiki/index.php/TRIPPEL_5pectrQgraph| 



0.0411" pixel ^ for both spectral regions while the spec- 
tral pixel size is 0.0105 A pixel-i and 0.0129 A pixel"! for 
the Ha- and Ca II IR - spectral regions respectively. To 
obtain a reasonable signal-to-noise ratio we observe with 
80 ms exposure times with both cameras. For each of 
the spectrogram cameras a slit-jaw camera was operated 
with the same exposure time and with a wideband filter 
centered close to the corresponding spectrogram wave- 
length. The filters used were centered at 6565 A and 
8714 A and had FWHM of 10 A and 100 A respectively. 
These slit-jaw images were used for the destretching of 
the spectrograms, see §[3l Furthermore, a slit-jaw cam- 
era was operated with an exposure time of 8 ms to obtain 
reference slit-jaw images less blurred by seening. The fil- 
ter used was centered at 6565 A with a FWHM of 10 A. 
To get context images for the spectrograms we obtain 
co-temporal and co-spatial observations of both the Ha 
line center and Ha continuum using the DOT. 

On 2006 May 04 we used a smaU pore in the NOAA 
active region 10878 (N14, E04, = 0.96) as adaptive 
optics lock point for the SST, see Fig.l, to obtain a time 
series of about 40 minutes (08:20:33-09:01:14UT) during 
good to excellent seeing conditions. During this time 
period we observe in "save all" mode, which in our case 
means a cadence of ~ 0.5 seconds. All data presented in 
this paper are from this time series. An overview of the 
observations can be seen in Fig[TJ 

3. DATA REDUCTION 

3.1. Aligning, destretching, and Fourier filtering 

Flat field and dark current images were constructed 
from the mean of 500 images, fol lowi ng the same proce- 
dure as described in (L angangen et ani2007[ ) . After these 
corrections, the spectrograms are aligned to the corre- 
sponding wide band slit-jaw images. This is done by 
correlating the continuum intensity in the spectrogram 
and the intensity close to the slit in the slit-jaw image. 
Furthermore, the continuum images from the DOT are 
aligned with the slit-jaw images, and since the Ha narrow 
band images on the DOT are aligned with the continuum 
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images we get a series of aligned spectrograms and both 
continuum and narrow band images, as seen in Fig[TJ 

Since the AO is not able to compensate for all the 
modes in the seeing, especially since the exposure times 
are 80 ms, post processing of the spectrograms as well 
as the images is desirable. Both the narrow-band and 
wide-band DOT images are po s t-pro cessed using the 
speckle method ()von der Luehd I1993D . The slit-jaw 
images obtained with 8 ms exposure times are post- 
processe d using the multi fra me blind deconvolution 
(MFBD, Ivan Noort et al.ll2005D . Due to the lack of spec- 
tral information outside the slit of the spectrograph, a 
similar reconstruction of the spectrograms is not possi- 
ble. A partial one dimensional compensation is, however, 
possible. The amount of destretching is obtained from 
the slit-jaw image with the same exposure time as the 
spectro gram. For th i s purp ose we use the routines devel- 
oped by I Shine et al.l (|1994f) . The amount of destretching 
from the tile covering the slit is then used to destretch 
the spectrogram. Due to the one dimensionality of this 
procedure it is clearly not a complete restoration method, 
hence the improvements are quite small. 

The spectrograms have a signal-to-noise ratio of typ- 
ically 180 at continuum level. To reduce the noise level 
we apply a conservative low pass Fourier filter technique. 
This reduces the noise significantly (about a factor two), 
without losing any significant signal. 

3.2. Wavelength calibration 

A mean solar spectrogram is constructed by adding 
the flat field spectra, altogether 500 exposures. The 
mean aberration co rrected spectrum is the n compared 
to the FTS atlas of lBrault fc Neckell (|1987f) . This atlas 
has proved to be well calibrated in wavelength and it 
shows no systematic offset in line shift s with wavelength 
(|Allende Prieto fc Garcia Loped [T9981) . Since the solar 
atlas is corrected for the Earth's rotation, the Earth's 
orbital motion and the Sun's rotation we automatically 
get a corrected spectrogram when we calibrate using the 
FTS atlas. Note that in this way our wavelength scale 
is the same as that of the FTS atlas; we therefore have 
to correct for the gravitational redshift and the convec- 
tive blueshift. The convective 'blueshift' is in fact a red- 
shift not caused b y classical convect ion in the case of 
the Ca II IR lines (|Uitenbroekll2006f ). We use the rest 
wave l ength given by the VALD database ([Pis kuno v et al.1 
fl995t iKupka et all [l99l lRvabchikova"etaIlll999l) . We 
correct for the gravitational redshift of 634 ms~^. The 
residual redshift in the mean profile of 425 ms^^ is 
caused by the intensity and velocity variations in the 
chromosphere, and poss ibly a lso by errors in the rest 
wavelength (.UitenbroeS 120061 ). No correction of this 
residual redshift is done. 

The slit covers a network region of about 38" or ap- 
proximately 28 Mm on the Sun. We have tested the 
difference between the "quiet" solar atlas and our "ac- 
tive region" mean spectrum . Even though the "active 
region" mean spectrum is quite different with respect to 
the shape of the line profile, the mean line center shift 
differs only with about 70 ms^^. Since time dependent 
variations in the wavelength calibration are known to be 
quite small we do not use any time dependent wavelength 
calibration. All in all we estimate the wavelength calibra- 
tions to have an accuracy of about 100 ms^^, the errors 



due to uncertainties in the rest wavelength are, however, 
of several hundred ms~^. 
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Fig. 2.— Velocity-time plot for the Ca II 8662 A line. Several 
DFs are seen as diagonal dark components across the spectral line. 



4. OBSERVATIONAL RESULTS 

Since the Ha line is a very wide spectral line with a flat 
center, it is not very sensitive to Doppler shifts. In addi- 
tion, the difficulties of understanding the line formation, 
make Ha less than ideal for measuring Doppler shifts and 
interpreting the results. Instead, we use the Ca II 8662 A 
line to measure Doppler shifts, only using the Ha line as a 
reference when we identify DFs. The reader should note 
that the Ca II 8662 A line also has several disadvantages. 
Firstly, it has a Fe blend in the blue wing, at 8661.9 A, 
that can interfere with large blueshifts in the line cen- 
ter. Secondly the line is very dynamic and shows several 
absorption and emission components, which complicates 
Doppler measurements. This line is nevertheless suited 
to measure the Doppler shifts in the DFs, both because 
it is quite narrow with steep wings and because DFs are 
easily identified as dark components moving across the 
line, as seen in Figl2l 

Since the spectrograms suffer from residual seeing ef- 
fects, the solar features (e.g. DFs, bright points etc.) wiU 
move somewhat both across and along the spectrograph 
slit. To reduce the effect of this kind of movement we 
re-bin the spectrograms to a spatial pixel-size of 0.12", 
and use only the highest contrast images in 6 seconds 
bins for this analysis. 

To visualize the DFs in a simple way and to make the 
identification of DFs easier, we make velocity-time (VT) 
plots, see Figl2] The DFs are identified by the diagonal 
dark components seen in the VT plots. To avoid miss- 
identifications of DFs we only identify DFs above the 
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Fig. 3. — Doppler velocities measured in 26 observed DFs 
(crosses). The least square linear fit to these points is plotted 
as a solid line in each panel. 

plage region in the lower part of the continuum images. 
The narrow band images show that this region clearly 
contains DFs, but the spatial resolution of the narrow- 
band images is not good enough to measure the motion of 
independent DFs. Despite the lack in spatial resolution 
we look for a drop of the intensity in the narrow-band 
images to accept a DF identification. Furthermore, we 
look for a co-spatial and co-temporal component in the 
Ha spectrograms. In this way we manually identify 26 
DFs. 

We measure the Doppler shift in these DFs by fitting 
a 4th order polynomial to the dark component, first ap- 
pearing in the blue wing, then tracking it through the 
line center, before it appears in the red wing. The iden- 
tification of the components is also done manually. The 
resulting Doppler measurements of the DFs arc shown in 
FigH 

A linear polynomial is fitted to each of the observed 
DFs. Tests with the reduced x^~inethod show that this 
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Fig. 4. — Scatter plots of different observables. The linear least 
square fit to the data is plotted as a dashed line. The Spearman 
correlation coefficients for each scatter plot are sh own in the lower 
right c orne r of each panel. The da ta obtained by I Hansteen et al.l 
112003) and lDe Pontieu et al.l | |200^ are shown as an inverse gray- 
scaled density image in the lower panel. 

fit is a good fit if the error is estimated to about 2 kms~^. 
An error of 2 kms"'^ is reasonable for these observations. 

Accepting the linear fits we can derive a mean value 
for the deceleration of 89 ± 25 ms~^, for the maxi- 
mum velocity of 11.3 ± 3.8 kms""'^, and life times of 
217 ± 39 sees, with errors given in one standard devi- 
ation. These values are all on the low side of earlier 
observations, the reason for this discrepancy will be dis- 
cussed in detail in fjni It should also be mentioned that 
the maximum velocity in the receding phase of the DF 
life, seems to have higher values than in the emerging 
phase, on average 2.4 kms""*^ higher. This result may be 
caused by the Fe blend in the blue wing of the Ca line. 

We get a good correlation between the maximum ve- 
locity and the deceleration, see FigU] lower panel. Using 
a Sp earman rank correlation test (see e.g. [Press et alj 
Il992f ) we get a correlation coefficient of 0.82. With the 
null hypothesis of no correlation we get a significance of 
4.1ct for falsification, which is a fairly strong statistical 
significance. This correlation has been shown to be a 
signature of the DF shock theory (IHansteen et aLll2006t 
iDe Pontieu et al.ll2007t iHeggland et al.ll2007t ). The weak 
negative correlation between the lifetime and the decel- 
eration, see Figd] upper left panel and the weak positive 
correlation between the maximum velocity and the life- 
time, seen in the upper right panel of Fig[¥]have a statis- 
tical weak significance of only 1.4cr and 0.9a respectively. 
Simila r weak correlations were found in lDe Pontieu et al.l 
((2001 . 

5. ERROR ANALYSIS 

The observations presented in this paper differ from 
previously presented observations of DFs in some funda- 
mental ways. Our data contain the full spectroscopic in- 
formation, while previous results have usually been based 
on narrow band images from one or a small number of 
spectral positions. This makes our observations more ro- 
bust when it comes to Doppler velocity measurements 
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compared to other methods. Unfortunately this advan- 
tage is only achieved at the expense of other disadvan- 
tages. In our case there is one main concern, namely 
the lack of spectral information from regions outside the 
slit. This lack of information leads to two problems in 
our analysis. Firstly parts of the DFs are missed, thus 
altering the correlations between the observables. Sec- 
ondly it makes us unable to post-process the observa- 
tions in a fully consistent way. Nevertheless these obser- 
vations contain important information of the DF dynam- 
ics. Since we can not correct for these shortcomings we 
will try to understand their effects on the observations. 

To get an estimate of the magnitude and effect on the 
observables due to the lack of information outside the 
slit, we investigate an Ha line- core imaging times-series , 
similar to the one presented bv lDe Pontieu et all (|2007D . 
We pick a random "slit" covering a similar plage region 
as in our spectroscopic observations. We identify 11 DFs 
crossing the slit, and measure the movement of these 
DFs. The part of the DFs outside the slit is removed to 
investigate the effect of the ID slit. Both the whole DF 
statistics and the reduced DF statistics show the same 
correlation between the deceleration and maximum ve- 
locity, but the correlation becomes weaker when parts of 
the DFs are removed. This indicates that the spectro- 
scopic results are still valid, but somewhat affected. The 
mean values are typically reduced with about 20%. We 
do expect this effect to be much less severe in Doppler 
measurements, since we are biased toward DFs that have 
strong LOS components, and hence are smaller in the im- 
age plane compared to the typical DF as seen in a narrow 
band image. Doppler measurements are usually strongly 
affected by the point spread function (PSF) since it con- 
taminates the resolved spectra with other spectra origi- 
nating from outside the position of the slit. In the present 
case the PSF is actually making the Doppler measure- 
ment more robust. The fact that the high Doppler shifts 
are seen as a separate component results in a small ef- 
fect from the PSF on these Doppler measurements. This 
is because at the wavelength of the separate component 
other spectra (from neighboring spatial pixels) usually 
are in the far wing. The low gradient of the far wing re- 
sults in a small effect on the Doppler component. Since 
this effect also works vice versa the high Doppler shifts 
will show in spectra where it otherwise would not, hence 
increasing the chances of observing the whole DF. 

6. NUMERICAL EXPERIMENTS 

To explain the lower velocities and decelerations in the 
present observations, realistic simulations arc needed. 
The simulation presented by Hanstccn ct al. (200^) re- 
produce some of the key results from their observations, 
and it should therefore be well suited for our purpose. 
We solve the radiative transfer equations for the atmo- 
spheres given by these simulations. This is done by using 
the MULTI code (Carlsson 1986). The simulated atmo- 
spheres are given in two dimensions, but the problem is 
simplified by solving for each column independently, each 
resolving 32.5 km of the solar surface. Each column is 
solved in full NLTE for the Ca line, while we include the 
line blend in LTE, but with a scattering term included 
in the source-function. When solving the radiative trans- 
fer CRD is assu med; the results w ill not be significantly 
ahered by PRD (jUitenbroeklll989f) . 
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Fig. 5. — Atmospheric parameters for three different times in 
the Hfe of a DF, 20 seconds (panel a), 100 seconds (panel b) and 
160 seconds (panel c) are seen. Velocity (solid line), temperature 
(dotted line), logarithm of the density (dashed line), and finally the 
contribution functions to relative intensity (dashed-dotted line) at 
the Doppler shifted wavelength are plotted. 

Since the simulations have several drawbacks, such as 
instantaneous hydrogen ionization, two dimensional ge- 
ometry, and lack of global 5-minute oscillations, it can 
not reproduce all details of the DFs. The simulations 
can nevertheless be used to illustrate the general physics 
involved. For this purpose we have chosen a DF-like fea- 
ture in the simulations, which is analyzed in detail. The 
DF-like feature is identified in a region in the simula- 
tions where we have magnetic fields emerging from the 
photosphere. In this region the fluxtubes are fairly ver- 
tical, which is similar to the expected inclinations in the 
observations. Several DF-like features can be identified 
in this region of the simulations, and most of them show 
si milar behavior a s the o ne we investigate in det ail. 

iHansteen et"all (|2006l ) and iDe Poiitieu et all (|2007l ) 
measure the movement of a bright top, which is caused 
by the steep temperature increase in the transition re- 
gion where it is believed that Ha has some contribution. 
Hence they measure the movement of a point in the im- 
age plane, that is well defined due to the big contrast seen 
between the top of the DF and the surroundings. This 
observational method is compared with the spectroscopic 
method. 

Since the calculation of the response function to ve- 
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locity (RFV) is based on linear approximations it is not 
possible to calculated RFV for DFs, which are inherently 
non-linear in nature. It is, however, possible to calculate 
the contributio n functions to relative absorption (CFRA) 
(jMagainl l 19861. The CFRA are fairly similar to the RFV 
since the Doppler shifts are seen as an absorption com- 
ponent. We calculate the contribution functions at the 
Doppler shifted wavelengths as measured in the different 
parts of the DF. The resulting contribution functions and 
the atmospheric velocity, temperature and density can be 
seen in FiglS] 

The wave seen in panel A in FiglS] steepens as the at- 
mospheric density drops about 4 magnitudes from the 
photosphere into the chromosphere. In the first stage of 
this DFs lifetime the transition region is found at about 2 
Mm. The Ca IR line samples a region below the maximal 
velocity, which is seen close to the transition region. At 
this time in the DF the measured Doppler shift is about 
6 kms^^, while the velocity close to the transition region 
is about 10 kms~^. This difference in velocities is very 
similar to the differences between the observations based 
on the motion of the TR and those based on the Doppler 
shifts, see Fig HI 

In panel B in Fig [5] the atmospheric parameters at a 
later time in this DFs lifetime are plotted. The transi- 
tion region has been pushed upwards by the shock, and 
it is now at a height of about 3.6 Mm. Even though the 
formation height of the Ca IR line is somewhat higher 
in the atmosphere at this stage of the DF lifetime, it is 
still formed below the transition region. The formation 
height of the Ca IR line is basically governed by the in- 
tegrated density above the formation height, hence the 
reason for the higher formation height at this stage of the 
DFs life time is the higher densities at equal geometric 
heights. The extension in height of the CFRA is of some 
significance to the measured Doppler shift, which in this 
case is close to zero. 

In panel C in Fig[5]the atmospheric parameters at the 
end of the DFs lifetime is plotted. The transition region 
height has declined, and is now seen at about 2.8 Mm. 
The atmospheric velocity close to the transition region is 
about —13 kms^^, while the measured Doppler velocity 
is about —8 kms^^. It should be noted that due to the 
lower height of formation in the Ca IR line the shocks 
sometimes overlap more than what is seen in the tran- 
sition region. This leads to sometimes shorter lifetimes 
in the DFs seen in the Doppler grams compared to the 
transition region. Furthermore some of the movement of 
the transition region is apparent motion due to sideways 
motion of the DFs across the slit. 

From the analysis above the intrinsic differences in the 
observational methods can explain the observed differ- 
ences in maximal velocities. The difference in decelera- 
tion can be explained in the shock framework, since the 
maximal velocities are lower, the deceleration must also 
be lower assuming similar lifetimes (as mentioned above 
the lifetimes are somewhat shorter in the spectra, but the 
relative reduction in lifetimes is smaller then the relative 
reduction in maximal velocities). 

The DFs in the simulated spectra have lower maxi- 
mal Doppler velocities, than in the observations. Typi- 
cal maximal values in the simulations are 9 kms~^ versus 
11 kms~^ in the observations. Furthermore the spread 
in the maximal velocities in the simulations are lower 



than in the observations. These differences are in part 
caused by a central reversal, which makes it impossible 
to measure the extreme Doppler shifts. The results are 
lower maximal velocities and life times, but no change 
in decelerations. These central reversals are caused by 
the simplified two dimensional geometry, which neglects 
the horizontal energy dissipation, hence increasing the 
temperature in the DFs. Tests from a three dimensional 
simulation show that this can account for some of the 
deficit in velocity, but not all. Other reasons can be dif- 
ferences in inclination angle of the local magnetic field in 
the observations and simulations. The shorter life times 
in the simulated DFs, due to lack of global 5-minutes os- 
cillations in the simulations will also affect the maximal 
velocities. 

7. CONCLUSIONS 

In this work we have presented co-spatial and co- 
temporal narrow band Ha images and spectroscopic mea- 
surements of the Ca II 8662 A line. These observations 
have been used to identify 26 DFs, and measure their 
Doppler shifts. A reduced-x^ analysis shows that the 
time evolution of the Doppler shifts are well approx- 
imated by a linear fit, if the measurement errors are 
about 2 kms~^. Using this approximation we derive 
values for the decelerations and maximal velocities for 
each DF. Scatter plots of the deceleration and maximal 
velocity show a strong positive correlation between the 
two. We also observe weak correlations between the de- 
celeration and lifetime and the maximal velocity and life- 
time. These results are supporting the shock- wave theory 
as explanation model for the DFs fHansteen et al.i r 2006l : 
iDe Pontieu~a l. 2007; Heggiand ct al. 2007 ). 

Furthermore, the Doppler shifts show that at least a 
subset of DFs are caused by mass moving up and down 
in the atmosphere. The values of the maximum velocity 
and decelerations are all somewhat lower than earlier re- 
ported values. Using numerical experiments we have ex- 
plained the differences in the two observational sets with 
the intrinsic differences in observational methods. Earlier 
observations have used the high contrast seen between 
the top of the DF and the background for measuring the 
proper motion of the DF. This high contrast is caused 
by the intensity increase due to contributions from the 
transition region. In the present observations the DF 
motion is measured using Doppler shifts, which are af- 
fected by the atmospheric conditions over the formation 
height. The formation height of the Ca II 8662 A line is 
much lower than the transition region. Since the shock 
amplitude is increasing from the formation height for the 
Ca IR line to the transition region we necessarily mea- 
sure lower velocities using spectroscopy. The difference 
in maximal velocities derived from the two methods in 
our simulations is about a factor two, which is about the 
same as observed. 

The one dimensional nature of the slit spectrograph 
somewhat affects our results, but experiments with nar- 
row band images show that this is not altering the results 
significantly. 
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7 



This research was supported by the European Commu- 
nity's Human Potential Program through the TOSTISP 
(contract HPRN-CT-2002-00310) program and by The 
Research Council of Norway through grant 170935/v30 
and through grants of computing time from the Pro- 
gramme for Supercomputing. The Swedish 1-m Solar 
Telescope is operated on the island of La Palma by the 
Institute for Solar Physics of the Royal Swedish Academy 



of Sciences in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto dc Astroffsica dc Canarias. 
The Dutch Open Telescope, operated by Utrecht Uni- 
versity at the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astroffsica de Canarias. 
This research has made use of NASA's Astrophysics Data 
System. 



REFERENCES 



Allcndc Pricto, C, & Garcia Lopez, R. J. 1998, A&AS, 129, 41 
Beckers, J. M. 1968, Sol. Phys., 3, 367 
Bel, N., & Leroy, B. 1977, A&A, 55, 239 

Brault, J., & Neckel, H. 1987, Hamburg Observatory anonymous 

ftp: ftp.hs.uni-hamburg.de 
Carlsson, M. 1986, A Computer Program for Solving Multi-Level 

Non-LTE Radiative Transfer Problems in Moving or Static 

Atmospheres (Uppsala Astronomical Observatory: Report No. 

33) 

Christopoulou, E. B., Georgakilas, A. A., & Koutchmy, S. 2001, 

Sol. Phys., 199, 61 
De Pontieu, B., Erdelyi, R., & De Moortel, I. 2005, ApJ, 624, L61 
De Pontieu, B., Erdelyi, R., & de Wijn, A. G. 2003, ApJ, 595, L63 
De Pontieu, B., Erdelyi, R., & James, S. P. 2004, Nature, 430, 536 
De Pontieu, B., Hansteen, V. H., Rouppe van der Voort, L., van 

Noort, M., & Carlsson, M. 2007, ApJ, 655, 624 
de Wijn, A. G., & dc Pontieu, B. 2006, A&A, 460, 309 
Hansteen, V., Dc Pontieu, B., Rouppe van der Voort, L., van Noort, 

M., & Carlsson, M. 2006, ApJ, 647, L73 
Heggland, L., Dc Pontieu, B., & Hansteen, V. H. 2007, ApJ, 666, 

1277 

Koza, J., Siitterlin, P., Kucera, A., & Rybak, J. 2007, in ASP Conf. 

Ser. 368, ed. P. Heinzcl, I. Dorotovic, & R. Rutten 
Kupka, F., Piskunov, N., Ryabcliikova, T. A., Stempcls, H. C, & 

Weiss, W. W. 1999, A&AS, 138, 119 
Langangcn, 0., Carlsson, M., & Rouppe van der Voort, L. 2007, 

ApJ, 655, 615 
Magain, P. 1986, A&A, 163, 135 
Michalitsanos, A. G. 1973, Sol. Phys., 30, 47 
Nishikawa, T. 1988, PASJ, 40, 613 

Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & 
Jcffcry, C. S. 1995, A&AS, 112, 525 

Press, W. H., Teukolsky, S. A., Vctterling, W. T., & Flanncry, 
B. P. 1992, Numerical recipes in FORTRAN. The art of scientific 
computing (Cambridge: University Press, — cl992, 2nd ed.) 



Rimmcle, T. R. 2000, in Proc. SPIE Vol. 4007, p. 218-231, Adaptive 

Optical Systems Technology, Peter L. Wizinowich; Ed., ed. P. L. 

Wizinowich, 218-231 
Rouppe van der Voort, L. H. M., De Pontieu, B., Hansteen, V. H., 

Carlsson, M., & van Noort, M. 2007, ApJ, 660, L169 
Rutten, R., Hammerschlag, R., Bettonvil, F., Siitterlin, P., & de 

Wijn, A. 2004, A&A, 413, 1183 
Ryabchikova, T., Piskunov, N., Stempels, H., Kupka, F., & Weiss, 

W. 1999, Physica Scripta, T83, 162 
Scharmer, G., Bjelksjo, K., Korhonen, T., Lindberg, B., & 

Petterson, B. 2003a, in Proceedings of the SPIE, Volume 4853. 

pp. 341-350 

Scharmer, G. B., Dettori, P. M., Lofdahl, M. G., & Shand, M. 

2003b, in Innovative Telescopes and Instrumentation for Solar 

Astrophysics. Edited by Stephen L. Keil, Sergey V. Avakyan . 

Proceedings of the SPIE, Volume 4853, pp. 370-380 
Shine, R. A., Title, A. M., Tarbell, T. D., Smith, K., Frank, Z. A., 

& Scharmer, G. 1994, ApJ, 430, 413 
Sucmatsu, Y. 1990, in LNP Vol. 367: Progress of Seismology of the 

Sun and Stars, ed. Y. Osaki & H. Shibahashi, 211 
Suematsu, Y., Wang, H., & Zirin, H. 1995, ApJ, 450, 411 
Tsiropoula, G., Alissandrakis, C. E., & Schmicdcr, B. 1994, A&A, 

290, 285 

Tziotziou, K., Tsiropoula, G., & Mein, P. 2004, A&A, 423, 1133 
Uitenbrock, H. 1989, A&A, 213, 360 
— . 2006, ApJ, 639, 516 

van Noort, M., Rouppe van der Voort, L., & Lofdahl, M. G. 2005, 

Sol. Phys., 228, 191 
van Noort, M. J., & Rouppe van der Voort, L. H. M. 2006, ApJ, 

648, L67 

von der Luehe, O. 1993, A&A, 268, 374 



